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ABSTRACT 

We discuss the emission and transport of polarized radio-band synchrotron 
radiation near the forward shocks of young shell-type supernova remnants, for 
which X-ray data indicate a strong amplification of turbulent magnetic field. 
Modeling the magnetic turbulence through the superposition of waves, we calcu- 
late the degree of polarization and the magnetic polarization direction which is 
at 90° to the conventional electric polarization direction. We find that isotropic 
strong turbulence will produce weakly polarized radio emission even in the ab- 
sence of internal Faraday rotation. If anisotropy is imposed on the magnetic-field 
structure, the degree of polarization can be significantly increased, provided in- 
ternal Faraday rotation is inefficient. Both for shock compression and a mixture 
with a homogeneous field, the increase in polarization degree goes along with a 
fairly precise alignment of the magnetic-polarization angle with the direction of 
the dominant magnetic- field component, implying tangential magnetic polariza- 
tion at the rims in the case of shock compression. We compare our model with 
high-resolution radio polarimetry data of Tycho's remnant. Using the absence 
of internal Faraday rotation we find a soft limit for the amplitude of magnetic 
turbulence, 5B < 200 /xG. The data are compatible with a turbulent magnetic 
field superimposed on a radial large-scale field of similar amplitude, 5B ^ Bq. 
An alternative viable scenario involves anisotropic turbulence with stronger am- 
plitudes in the radial direction, as was observed in recent MHD simulations of 
shocks propagating through a medium with significant density fluctuations. 



Subject headings: acceleration of particles, cosmic rays, methods: numerical, 
shock waves, supernova remnants, turbulence 
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Introduction 



The origin of Galactic cosmic rays and the mechanisms of their acceleration are among 
the most challenging problems in astrophysics. Shell-type supernova remnants (SNR) have 
long been thought to be the sources of cosmic rays, primarily through acceleration at their 
powerful forward shocks. Particle acceleration at coUisionless shocks is intrinsically effi- 
cient (e.g. Kang & Jones 2005) and arises from pitch-angle scattering in the plasma f lows 



that have systematically different velocities upstream and downstream of the shock (IBell 



19781 ). Detailed studies show that the acceleration efficiency and the resulting spectra de- 



of magnetic turbulence near the shock (e.g. 


Giacalone & JokiDii 


1996; 


2001; 


Giacalone 


2005; 


Bvkov &; Toptve:in 


2005h. part of which is self- 



am- 



plitude of the turbulence also sets the scale for the maximum energy to which a remnant 
may accelerate particles. For typical interstellar magnetic field values, SNRs can at best 
accele rate particles to 10^^ eV, w here the cosmic-ray spectrum shows a break known as the 
knee ( iLagage fc Cesarskvlll983al lbll If the cosmic rays would drive a turbu lent magnetic 
field to an amplitude m uch larger than the homogeneous interstellar field (ILucek fc Bell 



2000; iBell fc Luce 



( IVladimirov et al 



20061). 



20011) . particle acceleration may be faster and extend to higher energies 



Belli ( I2OO4J ) found that, rather than resonant Alfven waves, the current carried by drifting 
cosmic rays should efficiently excite non- resonant, nearly purely growing modes of magnetic 
turbulence on spatial scales much smaller than the cosmi c-ray Larn i or radius. MHD simula- 



tions that assume the co smic-ray current to be constant (IBelll 12004 . |2005| ; IZirakashvili et al. 



20081 ; iReville et al.ll2008l ) indeed indicate a strong magnetic-field amplification following an 
approximately isotropic plasma filamentation in the non-linear stage. On the other hand, 
recent kinetic simulations suggest that the amplitude of the field perturbations saturates 
at approximat ely the amplitude of the homogeneous upstream field on account of nonlinear 
backreactions (INiemiec et al.ll2008l ). Density fiuctuations in the upstream r egion will also 



disto r t the shock, leading to turbulent magnetic- field grow th downstream (IBalsara et al. 
2OO1I ; iGiacalone fc Jokipiill2007l ; Izirakashvih fc Ptuskinlljoosl ). 



While its micro-physics is not fully understood, magnetic-field amplification is observa- 
tionally required to explain that a large fraction of the non-t hermal X-ray emission on the 



rims of young SNRs is concentrated in narrow filam ents (e.g. iHughed |2000| ; iGotthelf et al. 



2OOII ; iHwang et al.l |2002| ; iBamba et al.l |2003| . |2005| ) . These filame nts can be interpreted 



either as lirn i ted b y rapid energy losses of the radiating electrons (jVink fc LamingI 12003 



Bamba et al.l 120031 ) or, alternatively, as c aused by rapid d amping of strong magnetic tur- 
bulence downstream of the SNR shock ( Pohl et "all I2OO5] ) . Both interpretations require 
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magnetic-field amplification, and the turbulently amplified field is expected to have a small 
wavelength. A very strong magnetic field is also suggested by time-variabil ity of patches of 



non-t hermal X-ray emission near the forward shock of SNR RX J1713-3946 (lUchivama et al. 



2007 ). although observational limits to the radio e niission of secondary elec trons (iHuang &: Pohl 
20081 ) and the cosmic-ray e/p ratio in general (IKatz fc Waxmaru 120081 ) indicate that we 



could observe jus t the build-up and decay of a magnetic structure (e.g. iButt et al.l 12008 
Bvkov et al.ll2008h . 



Because of its role in particle acceleration it is of paramount importance to understand 
the properties of strong magnetic turbulence near the forward shocks of SNRs. A signifi- 
cant uncertainty in our interpretations arises from the unknown spatial distribution of the 
strong, turbulent magnetic field. The properties of the turbulent magn etic field can be traced 
through synchrotron intensity fiuctuations (e.g. ICho fc LazarianI 120021 ) or through polarime- 
try. Here we present a search for signatures of such turbulence in polarized radio synchrotron 
emission. For that purpose we create parametrized models of magnetic turbulence and cal- 
culate the emission and transport of linearly polarized synchrotron radiation through the 
turbulent downstream medium, including the effects of Faraday rotation. The results are 
comp ared with pubhshed radio pola rimetry data of Ke pler's SNR (SN 1604) (IDeLaney et al. 
2Qm and Tycho's SNR (SN 1572) Jpickel et al.lll99lh . 



Since the degree of linear polarization and the polarization angle depend mainly on the 
structure of the turbulent magnetic field, our results are very generic and can be widely 
applied. The amplitude of the magnetic field enters only through a scaling parameter that 
describes the extent of depolarization through Faraday rotation. 



The model 



2.1. Parametrization of the turbulent magnetic field 



A turbulent magnetic field can be constructed via superposition of many magnetic waves 
with random orientation. A single magnetic wave carries a magnetic-field vector of the form 

B = Bosin(k-x) (1) 

and is assumed to fill all space. In order to ensure t hat V • B = 0, we choose rn agnetic-field 
components of the form (iGiacalone fc Jokipiilll994J : iNiemiec fc Ostrowskil l2006l ) 



By 



Bo,x sin {kyy + k^z + Cr,) 
Bo,y sin {kr,x + k^z + ay) 
5o,^ sin (k^x + kyy + a 



(2) 
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where the (7, terms denote randomly generated phase shifts. The projections of the wavevec- 
tor k are: 

kx = k sin 6 cos rj 

ky = ksinOsmrj (3) 
kz ~k cos d 

where the angles 9 and 77 are randomly selected to generate isotropic turbulence. The 
wavenumber k is randomly generated with uniform distribution in In /c, such that 27r x 10~^ < 
/c < 47r X 10~^ inverse cell units. This allows a range of wavelengths between 5 and 1000 
cell units, or 20 — 4000 periods along the line of sight. The amphtude of the magnetic wave 
follows 



BQ{k) = Bo{kmin) 



k 



-1 



2 



kmin ^ k ^ kuiaxi ('^) 



h ■ 

where we consider two possible values for the power-law index, q. For a Kolmogorov spectrum 
q = 'b/'i and for a flat spectrum q = 1. For a model with a Kolmogorov spectrum, the rms 
magnetic amplitude is 5B = ^ {B^, + By, + B"^,) ^ 10/xG, while a flat-spectrum model 
results in SB 20//G. The flat-spectrum models involve more small-scale fluctuations than 
do the Kolmogorov models. We do not know the true shape of the turbulence spectrum 
in SNRs, but comparing the results for the two magnetic-fleld models permits at least a 
qualitative estimate of the expected results for arbitrary spectra. 

The projections of the magnetic-fleld vector are determined in the same way as for k: 

5o,x = Bq sin C cos ^ 

Bo,y =5osinCsin^ (5) 
Bo,z =-BoCOsC, 

except that ( and ^ are new angles, still randomly generated such that — 1 < cosC < 1 and 
< ^ < 27r. It is straightforward to develop a turbulent magnetic fleld from the superposition 
of 1000 waves by summing each component as follows 

1000 

Bx,tot = ^ Box,i sin {ky^i y + k^^i z + a^^i) . (6) 

i=l 

While it is important to ensure that V ■ B = for each wave, the final magnetic field 
components must also fluctuate along all three coordinate axes, otherwise there would be no 
fleld reversals along the line of sight and hence no turbulent Faraday rotation. Therefore, we 
rotate the coordinate system by 45 degrees around each axis, thus changing the coordinates 
(x, y, z) to (x', y', z') via 




x' = — p -1 - V2 -1 + V2 V2 X (7) 
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The rows of the transformation matrix give the unit vectors of the primed coordinate system, 
ex', Cy/, and e^/, and the columns contain the unit vectors of the old coordinate system, e^, 
By, and 62;. The projections of the magnetic-field vector on the new coordinate axes are 

Bx' = ■ {Bx^tot, By^tot, B;^^tot) (8) 

and analogously for By! and B^' . Here Eq. |6] is written in term of the primed coordinates 

1000 

Bx,tot = ^ Box,i sin {ky^i Cy • x' + k^^i ■ x' + (Ja;,i). (9) 

i=l 

and correspondingly for By^tot and Bz^tot- It is straightforward to see that B^' now depends 
on x' etc. 

In the primed coordinate system, a three-dimensional grid is defined with size 20x20x20000 
cells or 50x50x20000 cells. The longest dimension (20000 cells) corresponds to the line of 
sight, while the other two dimensions are in the plane of the sky. 



2.2. Transport of polarized radio emission 

We are interested in modeling the observed polarized radio emission near the rims of the 
remnant, where any signatures of the possibly warped contact discontinuity can be ignored. 
Figure [1] shows a sketch of the geometry. Radio observations of SNRs suggest that the 
synchrotr on intensity, and therefore the spectrum of radiating electrons, follows a power law 



ju oc u ° (lGreenll200ll ). The emissivity also depends on the perpendicular component of the 



magnetic field in the emission region: 

J, = K {a) Bl+^ z/-" with Bl = Bl + (10) 



wher e K{a) is a constant depending on the density of relativistic electrons (IRybicki fc Lightman 



I979I ) . For an isotropic ensemble of electrons the synchrotron emissivity has a degree of linear 



polarization 

a + 0/6 



which is independent of frequency (iLe Rouxlll96ll ). For typical values ofa; = 0.5orQ; = 0.7 



the initial degree of linear polarization is 70% and 72%, respectively. 

Radiotelescopes typically measure the electric-field direction of a linearly polarized wave, 
which reflects the direction of acceleration of the radiating electron. The e/ectnc-polarization 
angle is then rotated by 90° to obtain the ma^neizc-polarization angle, which reflects the 
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Fig. 1. — Sketch of the geometry. We discuss the transport of polarized radiation for a hne 
of sight (LOS) near the rim of the SNR, right behind the forward shock (FS) but outside 
the contact discontinuity (CD). 
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orientation of the perpendicular component of the magnetic field at the location of the 
electron. It is therefore practical to directly calculate the magnetic-polarization angle, which 
can be compared with radio polarimetry data, for which the 90°-rotation has already been 
performed. 



It is useful to describe polarized emission using a complex polarized intensity (IBurn 



19661) 



P, = h exp(2zx) , (12) 

the phase of which is twice the magnetic-polarization angle, x- The factor 2 in the argument 
of the exponential accounts for the indistinguishability of polarization angles x ^"^^ X + ^• 
The complex emissivity is then similarly defined as 

= Pju exp {2%^)) , (13) 

where the initial magnetic-polarization angle, if), reflects the orientation of the magnetic field 
perpendicular to the line of sight through the relations 

Bx' Byl 

cosip = —— sini/j = (14) 
B± B± 

Polarized radio emission will suffer Faraday rotation upon passage through a magnetized 
plasma. The amount of rotation is proportional to the square of the wavelength, and the 
proportionality constant is usually defined as the Faraday depth of the source. The final 
polarization angle 

X(A^)=^ + 0A^ (15) 

where (j) is the Faraday depth calculated as a pathlength integral along the line of sight of the 
local electron plasma frequency, Up^e, the electron gyrofreq uency, a;^, and the inclination angle 



between the line of sight and the magnetic-field vector, ^ ( lBurnlll966l : iBrentjens fc de Bruyn 



20051 ). 



Also, Tie denotes the density of free electrons and B\os = B^' is the magnetic-field component 
along the line of sight (positive when pointing toward the observer). In the absence of 
absorption the transport of polarized radio emission is then given by line-of-sight integration. 

e^{s) exp[2i4>{s)X^]ds (17) 











Llos 

p j„{s) exp [2t (0(s)A^ + ^{s))] ds 
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where Llos is the length of the hne of sight, here 20,000 cells. We calculate the observed 
degree of polarization as 

n \Pu\ I / ds pMs) exp [2z {(P{s)\^ + ^{s))] \ 
lu J ds j^{s) 

and the observed polarization angle according to Eq. [121 

It is worth noting three particular aspects of our method. 

1. The degree of polarization does not depend on the rms amplitude of the perpendicular 
component of the magnetic field, only on the spatial fluctuations and the orientation. 
It also does not depend on the spectral index, if we ignore the minimal variation of the 
intrinsic degree of polarization, p (cf. eq. [TTl) . In the absence of Faraday rotation, i.e. 
if |0(s)|A^ <C TT, also the line-of-sight component of the magnetic field, Byos, cancels 
out of Eq. [iHl which is then solely dependent on the structure of the turbulent field. 
The frequency cancels as well, and therefore we will find that at high frequencies the 
degree of polarization is inevitably small, if the magnetic field is fully and isotropically 
turbulent on small spatial scales. A significant homogeneous field or anisotropy of the 
turbulence, e.g. through shock compression, will increase the polarization fraction. 

2. Faraday rotation will generally lead to additional depolarization, if |0(s)|A^ > vr. Our 
examples assume numerical values of = 1 cm~^, a line-of-sight length Llos = 10 pc, 
and a mean turbulent magnetic-field strength of 6B ^ 10/iG for the Kolmogorov 
models and 6B ^ 20yuG for the flat-spectrum models. Variations in one variable can 
be fully compensated by variation in another variable, in particular the wavelength 
A. Faraday rotation may simply become important at somewhat higher (or smaller) 
frequency than shown in our figures. Because of the quadratic wavelength dependence, 
relatively small changes in wavelength can compensate for relatively large changes in 
the other variables. 

3. We have described the radiation transport along a single line of sight. Since we simulate 
a grid of 20x20x20000 cells or 50x50x20000 cells, we can also account for beam depolar- 
ization for a beam size of 20x20 or 50x50 cells, respectively. This beamsize corresponds 
to a linear resolution of 0.1% or 0.25% of the line-of-sight length. Llos is typically 
a fair fraction of the SNR radius (see Fig. [T]), and hence our beam size corresponds 
to somewhat less than 0.1% of the angular radius of the SNR. Most published radio 
observations have larger beams, and thus we underestimate the beam depolarization. 
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3. Results 

3.1. Isotropic strong turbulence 

The most generic case to consider is isotropic turbulence without a large-scale homo- 
geneous magnetic-field component. The calculations should broadly describe situations in 
which 6B ^ Bq, as is often invoked for young SNRs. Figure [2] shows the degree of polar- 
ization as a function of frequency for a selection of magnetic-field models with Kolmogorov 
and fiat turbulence spectra. In all cases the polarized intensity was "beam- integrated" over 
an area of 20 by 20 cells. We find the behavior of polarized radio-synchrotron emission well 
characterized by the following statements: 

• Generally, magnetic-field models based on a Kolmogorov turbulence spectrum tend to 
give a higher polarization degree than models for fiat turbulence spectra. This behavior 
can be understood with a toy model, in which the magnetic field is assumed constant in 
a zone of "wavelength" /, but having random variations from zone to zone. If there are 

zones on the line of sight, the mean observed degree of polarization would be about 
70% /^/N. For the Kolmogorov models we find an observed degree of polarization 
around 8%, or N = 80, or / ~ 250 cells. For fiat-spectrum turbulence, which has 
much higher amplitudes at small scales, we find 4%, or N = 300, or / ~ 65 cells. We 
can extrapolate to turbulence of smaller wavelength than simulated here, for which we 
expect the observed degree of polarization to be less than 4% on average. 

• The observed polarization angles are uniformly distributed on scales larger than the 
characteristic wavelength of the turbulence, which is easily understood in the frame- 
work of the toy model described under the first bullet above. Because our turbulence 
model is built through the superposition of waves, two lines of sight close to each other 
will pass through magnetic-field structures that are somewhat correlated, whereas lines 
of sight far from each other will be fully uncorrelated. 

• Figure [3] shows the small-scale distribution of the observed magnetic-polarization angle 
at 4.8 GHz for a Kolmogorov model. To be noted from the figure is that the polarization 
angle varies by less than 0.3 radian over an area of 50x50 cells, whereas it varies 
wildly over significantly larger areas. Therefore, the beam depolarization is small, 
if we beam-integrate the polarized intensity only over the 50x50 cells, or 20x20 cells 
as done for Fig. [2l Real radio measurements have a spatial resolution that is much 
worse than 0.25% of the line-of-sight length (50 cells vs. 20,000 cells). Therefore, 
beam depolarization may be much more significant in the real data, and the observed 
polarization degree may be lower than shown in Fig. [2j 
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Fig. 2. — The degree of polarization as a function of frequency for a selection of six turbu- 
lence models assuming Kolmogorov spectra (top panel) and six models based on flat spectra 
(bottom panel). 
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0.485 0.512 0.543 0.575 0.609 0.645 0.683 

Fig. 3. — The spatial variation of the observed magnetic polarization angle at 4.8 GHz, where 
internal Faraday rotation is negligible. A turbulence model with Kolmogorov spectrum was 
chosen. The plot covers a surface area of 50x50 cells, while the line of sight is 20,000 cells 
long. 
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At higher frequencies the observed polarization degree is determined only by depolar- 
ization through the intrinsic magnetic-polarization angle, ip, which reflects the orienta- 
tion of the perpendicular magnetic-fleld component along the line-of-sight (cf. Eg. [T^ . 

In Fig. [2] we see substantial fluctuations in the degree of polarization between different 
magnetic-fleld models with the same turbulence spectrum. We expect that for a wider 
beam much of this variation would be smoothed out. 

Below about 1 GHz Faraday rotation becomes important. If the product of rms 
magnetic-fleld strength, 6B, free-electron density, rze, and line-of-sight length, Llos; in 
any SNR is different from the canonical values described in Sec. 12.21 then the onset of 
Faraday depolarization may occur at a somewhat lower or higher frequency. 

Because the magnetic fleld is turbulent, the eff ective Faraday rotation in th e observed 



polarization angle no longer follows a A^-law (IBrentjens fc de Bruyrul2005l ). Faraday 
rotation in the foreground - not considered here - will follow the A^-dependency, but it 
generally cannot be extracted using a A^-law because it is superimposed on the internal 
Faraday rotation. 

The strong Faraday rotation at low frequencies reduces the polarization degree only 
weakly, because the polarized emissivities are already fully randomized at high fre- 
quency. 



3.2. Mixtures of turbulent and homogeneous fields 

The magnetic fleld near the rim of SNRs may not be fully turbulent, for example if the 
magnetic-fleld ampliflcation saturates at an amplitude 6B ^ Bq. We can model a mixture of 
turbulent and homogeneous flelds by adding to our models of fully turbulent magnetic fleld 
a homogeneous fleld that is oriented in the x-direction. 

In Fig. m we show the degree of polarization as a function of frequency for mixtures 
of turbulent and homogeneous flelds of varying strength. The polarized intensity has been 
"beam-averaged" over 50x50 cells, or 0.25% of the line-of-sight length, and therefore the 
curve for a fully turbulent fleld is smoother at low frequency compared with those shown in 
Fig. [21 

Largely independent of the choice of turbulence model, a homogeneous magnetic fleld of 
the same strength as the rms turbulence amplitude, Bq = 6B, drives the degree of polariza- 
tion into the range 30% to 35%. A homogeneous magnetic fleld at i?o = 2 6B gives a degree 
of polarization at the level 55% to 60%. However, at low frequencies where internal Faraday 



0,6 F 




1 10 
Frequency (GHz) 



Fig. 4. — The observed polarization degree as a function of frequency for models with a 
homogeneous field component of varying strength (50%, 100%, or 200% of the rms turbulence 
amplitude, SB). The top panel shows results for Kolmogorov turbulence, and the lower panel 
displays the degree of polarization in the case of turbulence with flat spectrum. 
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Model 


Turb. only 


+50% 


+ 100% 


+200% 


Flat_01 


40.59 


8.42 


4.21 


2.58 


Flat_02 


0.43 


1.92 


1.20 


0.81 


Klm_01 


35.59 


22.52 


5.67 


1.85 


Klm_02 


87.57 


12.84 


4.13 


2.43 



Table 1: Magnetic polarization angles, given in degrees relative to the direction of the homo- 
geneous field (i?oex), for varying ratios Bq/{5B) and a selection of turbulence models. The 
frequency is 4.8 GHz, and so internal Faraday rotation is very small. 



rotation is important the degree of polarization remains below 10%, and for fiat-spectrum 
models, i.e. a very small effective turbulence wavelength, the polarization degree falls below 
3%. 

Table [T] gives the observed magnetic-polarization angles at high frequency for different 
relative amplitudes of the homogeneous field, which is oriented in the x-direction (zero angle). 
Recall that the magnetic-polarization angle differs by 90° from the electric polarization 
angle. Already for Bq = 6B the magnetic-polarization angles refiect the direction of the 
homogeneous field component to within a few degrees. At low frequencies, where internal 
Faraday rotation is strong, two lines of sight through statistically independent turbulence 
give a different polarization angle as if randomly selected, irrespective of the magnitude of 
the homogeneous field, but the degree of polarization is small. 



3.3. Shock-compressed turbulence 

If strong magnetic tu rbulence is generated upstream of the shock, for example by the 



streaming of cosmic rays (iBelll |2004| ). then it will change its properties when transmitted 
through the shock. Here we investigate this situation for the simplified case that the mod- 
ification by the shock can be fully described by a compression of the spatial scales in the 
direction of the shock normal (here the x-axis) and an increase of the perpendicular magnetic- 
field components. The change in spatial scale is achieved by reallocating the gridpoints of 
the mesh on which the magnetic field is calculated (see sec. 12. ip . Instead of sampling with 
Ax = Ay = Az the magnetic-field model is interpreted as sampled with Ax = {Ay)/n 
where k is the compression ratio. The compression of the perpendicular field components is 
effected by the substitution By — > nBy and likewise for B^. 

In Fig. [5] we show the polarization degree as a function of frequency for shock compres- 
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Model 


Turb. only 


K = 2 


K, = 3 


Flat_01 


40.59 


91.42 


93.48 


Flat_02 


0.43 


91.84 


92.63 


Klm_01 


35.59 


85.66 


89.46 


Klm_02 


87.57 


92.48 


93.72 



Table 2: Magnetic polarization angles, given in degrees relative to the shock normal, for 
varying compression ratios and a selection of turbulence models. The frequency is 4.8 GHz, 
and so internal Faraday rotation is very small. 



sion with k = 2 and k = 3, where the "beam-averaging" is effectively over 20x20 and 15x15 
cells, respectively, so the "beamsize" decreases with increasing k. Cosmic rays can modify 
the shocks at which they are accelerated, resulting in a reduction of the compres sion ratio 
at th e subshock of the thermal gas and an increase in the overall compression (IReynolds 



20081 ) . The choice of compression ratios in figure is solely meant to illustrate the behavior 



of the polarization, and in fact the changes in the degree of polarization will be only more 
pronounced for larger compression ratios. To be noted from the figure is that shock com- 
pression will substantially increase the degree of polarization at higher frequencies: already 
for a compression ratio k = 2 the degree of polarization falls into the range 35% to 40%, 
largely independent of the choice of turbulence model. As in the case of a mixture with a 
homogeneous field, Faraday rotation will strongly depolarize the emission at low frequency. 
The onset of efficient internal Faraday rotation is observed at somewhat higher frequency 
than in the case of a homogeneous field component perpendicular to the line of sight, be- 
cause the shock compression also increases the line-of-sight component of the magnetic field 
and hence the Faraday depth. While shock compression or a homogeneous field compo- 
nent can increase the degree of polarization, they will inevitably also modify the observed 
magnetic-polarization angle which is shown in Table [2l Even for a moderate compression 
ratio K = 2, at high radio frequencies the magnetic-polarization angles lie around 90° to the 
shock normal with very little scatter (the electric polarization is along the shock normal). 
At low frequencies Faraday rotation randomizes the polarization angle, but the degree of 
polarization is generally small. We can conclude that shock compression will not give a high 
degree of polarization > 20% and a magnetic polarization aligned with the shock normal. 
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Fig. 5. — The observed polarization degree as a function of frequency for models with shock- 
compressed turbulence for a choice of compression ratios, k. The top panel shows results for 
Kolmogorov turbulence, and the lower panel displays the degree of polarization in the case 
of turbulence with flat spectrum. 
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4. Discussion and summary 



We have seen that isotropic strong turbulence will produce weakly polarized radio emis- 
sion even in the absence of internal Faraday rotation. The polarization angle can be expected 
to vary on spatial scales of the order of the typical wavelength of the magnetic-field fluctu- 
ations, and therefore radio polarimetry data of very high angular resolution are needed to 
observe those variations. 

If anisotropy is imposed on the magnetic-field structure, the degree of polarization can 
be significantly increased, provided internal Faraday rotation is inefficient. Both for shock 
compression and a mixture with a homogeneous field the increase in polarization degree goes 
along with an alignment of the observed magnetic-polarization angle with the direction of 
the dominant magnetic-field components. In the case of shock compression we therefore 
expect tangential magnetic polarization at the rims of SNRs, or electric polarization vectors 
that are predominantly radial in orientation. 

Few young SNRs are suitable for the study of strong magnetic turbulence through ra- 
dio po larimetry. Some r emnants like Cas A feature ejecta beyond the nominal forward 



shock (IHwang et al.ll2004l ). and therefore our assumption of a simple, spherically symmetric 
structure does not apply. Also, the line-of-sight should not pass through the contact dis- 



continuity, at which th e magnetic field may be very strong (IRosenau fc Frankenthall Il976 



Lvutikov fc Pohibooi). It is well known that the contact discontinuity is hydrodynamically 



unstable (e.g. iBlondin fc EUisonl I2OOII ). and so in projection it will appear as an extended 



feature with a turbulent field structure. Also, if the SNR efficiently accelerates cosmic rays, 
the contact discontinuity will be cl oser to the forward s hock than in a purely hydrodynamical 
SNR. Using X-ray measurements, IWarren et al.l (120051 ) on average find traces of the contact 
discontinuity in Tycho's SN R out to 93% of the projected radius of the forward shock. 
Cassam-Chenai' et al.l (120081 ) use the same technique on data of the remnant of SN 1006 
and find that in the regions of bright non-thermal X-ray emission the contact discontinuity 
extends all the way to the forward shock. In both remnants the proximity of forward shock 
and the contact discontinuity presumably arises from a combination of hydrodynamical in- 
stabilities of the contact discontinuity and structural modifications on account of cosmic-ray 
acceleration. In any case, in SN 1006 we may not find a line of sight that is clearly inside 
the forward shock but outside the contact discontinuity. For Tycho's SNR radio-polarimetry 
data are required with an angular resolution around 1% of the angular radius, so the forward 
shock and the contact discontinuity are at least a few beamwidth apart. 



DeLaney et al.l (120021 ) observed the polarized synchrotron emission from Kepler's SNR 



at 6 cm and 20 cm wavelength. After rotating the measured electric polarization by 90°, 
they found predominantly radial magnetic polarization in the outer regions of the remnant. 
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where the degree of polarization was a few per cent at 6 cm and less than that at 20 cm. 
The angular resolution is moderate, though, and the beam size is about 7% of the projected 
radius of the forward shock. This is a factor of 100 worse than the beamsize assumed for 
the figures in this paper, and a direct comparison is therefore difficult. Also, the contact 
discontinuity and the forward shock are not clearly separated at this resolution, and so any 
inferred field orientation and turbulence level cannot be unambiguously associated with the 
magnetic-field structure directly behind the forward shock. 



Dickel et al.l (1l99ll ) presented radio polarimetry data of Tycho's SNR at 6 cm and 20 cm 
wavelength. The spatial resolution is about 0.7% of the forward-shock radius, roughly com- 
parable to 300 cells in our magnetic-field model. The radio morphol ogy can be well describe d 



by an outer rim and an inner shell of outer radius 0.92 SNR radii ( iKatz-Stone et al.ll2000l ). 
The outer rim is positionally coincident with the outer periphery of the X-ray emission, sug- 
gesting that the inner shell marks the location of the contact discontinuity and the reverse 
shock. We conclude that our model is applicable to polarized radio emission from the outer 



rim. 



In the outer rims, the percentage polarization at 6 cm wavelengt h is typically in th e 



range 20% to 30%, and the field orientation is radial (cf. Fig. 5a of iDickel et al.l Il99ll ). 
Beginning about 10" (or 7% of the SNR radius) toward the interior of the remnant, the degree 
of polarization is lower and the polarization angle varies on scales of about 10", consistent 
with large scale turbulence near the contact discontinuity. Internal Faraday rotation is 
likely negligible, because no deviation from a A^-law is observed between 6 cm and 20 cm. 
The lack of efficient internal Faraday rotation can be used to set limits on the amplitude 
of the turbulent magnetic field. The post-shock gas density in Tycho is not well known, 
but the estimates include ou r canonical numbe r Ug = 1 cm^^. The distance is likely in 



the range 1.5 kpc to 3 kpc (jSmith et al.lll99ll ). and therefore the line-of-sight length at 



the rim is Llos ~ 0-5 pc. For the canonical numbers used in our model internal Faraday 
rotation becomes important at about 1 GHz, which with the gas density and physical size 
of T ycho's SNR r e quires 6B ~ 200 fiG. If the field strength were significantly larger than 



that, iDickel et al.l (Il99ll ) should have observed internal Faraday rotation. 



Comparing with our model results presented in Sec. [31 we find no evidence for isotropic 
or shock- compressed magnetic turbulence that is amplified to an amplitude larger than the 
homogeneous magnetic field. The radio data are compatible with two scenarios: a turbulent 
magnetic field superimposed on a radial large-scale field of similar amplitude, i.e. a parallel 
forward shock with 6B ~ Bq, or, alternatively, anisotropic turbulence that has a somewhat 
larger magnetic-field amplitude in the radial direction. 

Anisotropic magnetic turbulence with preferentially radial orientation is not expected 
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in models involving magnetic-field amplification by streaming cosmic rays in the upstream 
region, in particular on account of shock compression. We note with interest that recent 
MHD simulations of magnetic-field growth following shock distortions on account of den- 
sity fiuctuations in the upstream medium suggest that the radial magnetic-field coni ponents 
may be slightly stronger than those in the shock plane ( Zirakashvili fc PtuskinI [20081 ). Addi- 
tional simulations with higher resolution are needed to confirm these findings and establish 
their dependenc e on the wavelength of the density fiuctuations and other parameters (e.g. 



Cho et all 120081 1. 
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